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►  Modeling  of  JP-8  fuel  reformer  integrated  with  a  tubular  solid  oxide  fuel  cell  (SOFC)  stack. 

►  Autothermal  reformer  (ATR)  is  characterized  to  obtain  optimum  performance. 

►  System  does  not  require  external  energy  and  steam  at  steady  state. 

►  System  efficiencies  of  39.5%  and  32.6%  observed  at  ATR  temperatures  of  700  °C  and  850  °C,  respectively. 

►  Optimum  operating  range  suggested  by  conducting  parametric  analyses. 
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An  onboard  autothermal  reformer  (ATR)  integrated  with  a  SOFC  stack  offers  potential  for  high  energy 
efficiency  and  utilization,  low  emission  and  quiet  operation  avoiding  cost  associated  with  hydrogen 
storage  and  infrastructure.  Such  a  system  can  be  a  viable  and  attractive  option  especially  for  military's 
need  for  quiet  and  less  pollutant  Auxiliary  Power  Unit  (APU)  and  Mobile  Electric  Power  (MEP)  units  in 
temporary  and  permanent  base  camps  [1,2],  A  5  kWe  Solid  Oxide  Fuel  Cell  (SOFC)  system  operating  on 
desulfurized  JP-8  fuel  was  modeled  using  Aspen  Plus  process  simulation  software  to  examine  the  effects 
of  oxygen  to  carbon  ratio  (O2/C)  and  steam  to  carbon  ratio  (H2O/C)  at  different  ATR  operating  temper¬ 
atures  (700—850  °C),  while  keeping  the  SOFC  stack  temperature  constant  at  910  °C.  Anode  recycle  steam 
and  heat  have  been  used  to  reform  the  desulfurized  JP-8  fuel  which  would  make  the  system  lighter  and 
compact  for  mobile  application.  The  system  modeling  revealed  a  maximum  net  AC  efficiency  of  39.5%  at 
700  °C  and  a  minimum  of  32.6%  at  850  °C  ATR  operating  temperatures,  respectively.  Sensitivity  analysis 
with  respect  to  fuel  utilization  factor  ( Uf )  and  current  density  (/)  were  also  conducted  to  identify  the 
optimum  operating  window. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFC)  are  high  efficiency,  low  emission, 
solid  state,  and  quiet  energy  conversion  devices.  But  commercial 
applications  are  still  limited  due  to  the  availability  of  fuel  such  as 
hydrogen  (H2)  and  high  initial  cost.  Integrating  a  hydrocarbon  fuel 
reformer  with  the  SOFC  system  can  minimize  the  costs  and  prob¬ 
lems  associated  with  hydrogen  storage  and  delivery.  The  United 
States  Army  is  interested  in  developing  a  fuel  cell  based  portable 
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power  system  to  support  their  energy  intensive  devices.  “Silent 
Watch”,  Army’s  term  for  very  deployable  power  systems  with  low 
acoustic  and  stringent  infrared  signature  levels,  make  it  impossible 
to  use  the  conventional  diesel  engines  due  to  the  loud  noise  in 
some  critical  operations  [1].  Moreover,  the  large  array  of  commu¬ 
nication  and  situational  awareness  electronics  used  in  fields  cannot 
be  supported  by  only  batteries,  for  long  durations,  making  this  need 
for  power  more  acute  [1],  These  types  of  systems  also  have  high 
potential  to  provide  power  in  the  temporary  and  permanent  mili¬ 
tary  base  camps  as  a  MEP  source  to  provide  power  for  mobile 
kitchen  units,  laundry  units,  heating  and  refrigeration  units,  well 
kit,  printing  plant,  hospital  maintenance,  bakery  plant,  water 
purification,  in  addition  to  support  the  combat,  communication  and 
missile  systems  [2],  As  JP-8/jet  fuel  is  the  main  logistical  fuel  of  the 
U.S.  Army,  an  APU  combined  with  JP-8  fuel  processor  can  provide 
a  viable  power  option.  It  can  also  be  modified  and  run  on  diesel  fuel 
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Nomenclature 

w 

component  width,  m 

yi 

molar  fraction  of  gaseous  component  i 

Roman  letters 

Greek  letters 

A,  B 

constants  in  Eqs.  (17)  and  (18) 

Af/ 

molar  Gibbs  Free  energy  of  formation,  J  mol"1 

D 

diameter  of  SOFC  cell,  m 

8o2 

constants  in  Eqs.  (24)  and  (25) 

Dik 

binary  diffusion  coefficient  of  species  i  and  k,  m2  s"1 

e 

electrode  porosity 

Dki 

Knudsen  diffusion  coefficient,  m2  s"1 

rizct 

activation  polarization,  V 

A(eff) 

effective  diffusion  coefficient  of  species  I,  m2  s"1 

Vcon 

concentration  polarization,  V 

E 

activation  energy,  J  mol-1 

’Jel.gross 

gross  AC  electric  power,  W 

F 

Faraday  constant,  C  mol-1 

57el.net 

net  AC  electric  power,  W 

I 

current,  A 

5?ohm 

ohmic  polarization,  V 

j 

current  density,  A  m  2 

V 

fuller  diffusion  volume 

K 

Pre-exponential  factor,  A  m“2 

? 

electrode  tortuosity 

l 

component  thickness,  m 

P 

resistivity,  Q  m 

LHV 

lower  heating  value  of  fuel,  J  mol-1 

Mi 

molecular  weight  of  gaseous  components  i,  kg  kmol"1 

Subscripts 

h 

molar  flow  rate,  mol  min-1 

a 

anode 

"fuel 

molar  flow  rate  of  fuel,  mol  min-1 

ac 

activation 

P 

pressure,  bar 

c 

cathode 

Pi 

partial  pressure  of  gaseous  components  i,  bar 

con 

concentration 

PsOFC 

SOFC  operating  pressure,  Pa 

e 

electrolyte 

Pel.DC 

DC  power,  W 

eff 

effective 

Pcl /vc 

AC  power,  W 

eq 

equivalent 

Pcomp 

compressor  power,  W 

i 

gaseous  component 

rp 

mean  radius  of  electrode  pore,  m 

int 

interconnector 

Rg 

universal  gas  constant,  8.314  J  mol"1  K 

k 

second  gaseous  component 

Pact 

specific  resistance,  Q  m2 

m 

mean 

Tavg 

average  temperature,  I< 

ohm 

ohmic 

Top 

SOFC  operating  temperature,  K 

Uf 

fuel  utilization 

Superscripts 

V 

voltage,  V 

a 

anode 

Vn 

Nernst  voltage,  V 

c 

cathode 

and  eventually  integrate  in  heavy  duty  trucks,  refrigeration  units, 
recreational  vehicles  (RVs)  due  to  their  long  inefficient  idling 
operations  [3],  It  can  also  be  an  attractive  option  for  aviation  APUs 
without  major  modification  of  the  existing  fuel  system  because  of 
lower  emission  and  usability  of  the  exhaust  heat  to  lower  the 
maximum  onboard  electric  load  [4], 

A  number  of  investigations  of  SOFC  based  processes  modeled 
using  Aspen  Plus  has  been  found  in  literature  using  different 
fuels.  Zhang  et  al.  [5]  used  natural  gas  in  their  study.  Doherty 
et  al.  [6]  used  two  types  of  biogas  to  investigate  system  perfor¬ 
mance.  Biomass-based  MWe  scale  combined  systems  have  also 
been  studied  comprising  of  gas  turbine  and  SOFC  module  to 
investigate  the  cycle  efficiency  and  optimum  operating  range  for 
stationary  applications  [7—10].  Baratto  et  al.  [11]  studied  a  diesel 
based  APU  for  commercial  heavy  duty  trucks  without  recycling  of 
the  anode  exhaust  heat  and  steam  and  showed  experimental  and 
modeling  agreement  of  results.  Roychoudhury  et  al.  [12]  investi¬ 
gated  a  5  kWt  JP-8  logistic  fuel  processor  with  external  steam 
generator.  They  mainly  focused  on  the  fuel  processor  and  did  not 
include  the  SOFC  stack  to  understand  the  overall  performance  of 
the  system.  Therefore,  a  complete  system  using  a  JP-8  fuel  based 
SOFC  with  APU,  which  would  address  the  steam  management 
needs  and  make  the  best  use  of  exhaust  heat,  is  not  available  in 
the  literature. 

The  objective  of  this  study  is  to  evaluate  a  5  kWe  tubular  SOFC 
based  system  that  operates  on  autothermally  reformed  desulfur¬ 
ized  JP-8  fuel  using  Aspen  Plus  process  simulation  software.  The 
adiabatic  ATR  was  simulated  and  characterized  from  700  °C  to 
850  °C  by  varying  the  H2O/C  and  O2/C  ratios  to  obtain  optimum  H2 


and  CO  contents  in  a  non-coking  system.  A  portion  of  anode 
exhaust  has  been  mixed  with  the  fresh  fuel  to  provide  the  steam 
and  heat  required  for  reforming.  Finally,  parametric  analyses  were 
conducted  by  varying  fuel  utilization  factor  (U/)  and  current 
density  (j)  to  identify  the  system’s  performance  at  different  load 
conditions. 

2.  Model  description 

2.1.  Model  assumptions 

The  model  is  based  on  zero  dimensional  approach  with  steady 
and  isothermal  operation.  All  working  fluids  are  considered  as  ideal 
gases  neglecting  any  pressure  drop.  It  is  assumed  that  reforming 
and  shifting  reactions  reach  chemical  equilibrium.  Only  H2  is 
assumed  to  be  electrochemically  oxidizable,  CO  is  assumed  to  take 
part  in  the  water  gas  shift  (WGS)  reaction  and  the  lower  hydro¬ 
carbons  (CH4,  C2H6,  C3H8  etc)  are  reformed  to  H2  [5,6], 

2.2.  Process  flow  model 

A  schematic  view  of  the  proposed  5  kWe  SOFC  system  is 
depicted  in  Fig.  1.  The  stream  named  ‘FUEL’  represents  the  fresh  JP-8 
fuel,  which  is  heated  to  its  higher  boiling  temperature  and 
compressed  in  the  compressor,  ‘CM-1’  block.  The  discharge  pres¬ 
sure  of  the  compressor  is  set  to  the  pressure  ratio:  Pfuei/PsoFC  =  3.0 
to  drive  the  recirculated  anode  gas  along  [5,6,13],  A  mixer  block 
named  ‘EJECTOR’  mixes  the  recirculated  anode  gas  with  the  fresh 
fuel.  A  particular  H2O/C  ratio  is  manipulated  by  changing  the  split 
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Fig.  1.  Process  flow  diagram  of  the  5  kWe  SOFC  system. 


fraction  using  a  Fsplitter  module  named  as,  ‘SPUTTER’.  H2O/C  ratio 
is  defined  as  the  molar  ratio  of  the  number  of  H2O  molecules  and 
the  combustible  C-atoms  entering  into  the  ‘REFORMER’  [13].  H2/C 
and  CO/C  ratios  are  defined  as  the  ratio  of  molar  flow  rate  of  H2  or 
CO  at  ATR  outlet  and  inlet  streams,  respectively.  A  design  specifi¬ 
cation  block  is  used  to  maintain  the  desired  H2O/C  ratio  by 
manipulating  the  split  fraction. 

A  Gibbs  reactor  module  named  ‘REFORMER’  is  used  to  simulate 
adiabatic  autothermal  reforming  of  JP-8  fuel.  A  heater  module 
named  ‘HEATR-1’,  preheats  the  air  coming  into  the  reformer  to 
supply  O2  flow  rate  in  the  reformer  to  make  it  adiabatic,  to  achieve 
the  maximum  reformer  efficiency  [  14],  The  governing  equations  of 
the  autothermal  reaction,  water  gas  shift  (WGS)  reaction  and 
steam  reforming  (SR)  reactions  are  shown  in  Eqs.  (l)-(3), 
respectively. 

C„Hm  +x(02  +3.76  N2)  +  (2 n  -  2x)H20^nC02 


+  (2n  -  2x  +  H2  +  3.75xN2 

(D 

co  +  h2o«co2  +  h2 

(2) 

CxHj,  +  xH20<-»xC0  +  (x  +  y/2)  H2 

(3) 

The  reformed  fuel  in  the  ATR  outlet  mainly  contains  H2  and  CO 
gases.  Lower  hydrocarbons  e.g.  CH4,  C2Fl6,  C3H8  can  also  be  present 
in  the  reformed  gas  and  assumed  to  take  part  in  the  internal 
reforming  reactions  (Eq.  (3))  in  the  ‘ANODE’  chamber.  The  elec¬ 
trochemical  reactions  are  as  follows: 

Cathode  half  -  cell  reaction  :  1/2  02  +  2e  ->02  (4) 


named  ‘CM-2’  and  later  heated  in  the  heater  block  ‘HEATR-2’.  The 
pre-heated  air  is  further  heated  in  the  heat  exchanger  module 
(HeatX)  named,  ‘RECPT’  by  utilizing  the  hot  exhaust  gas  stream 
coming  out  from  the  post  combustor,  ‘POST-CMB’,  which  simulated 
by  a  stoichiometric  reactor  and  ‘HEATR-5’.  The  compressed  and 
heated  air  (stream  ‘14’)  goes  into  the  ‘CATHODE’  block,  which  is 
a  separator  module  that  separates  02  from  air  needed  by  the 
electrochemical  reaction  and  consequently  simulates  the  cathode 
of  SOFC.  Aspen  Plus  calculates  the  flow  rate  of  02  based  on  the 
equivalent  H2  flow  rate  at  the  anode  as  follow: 

rcH^eq  =  hH2  ATR  +  l(nCOATR)  +  4(fiCH4  ATR)  +  7(iiC2H6  ATR) 

+  10(nC3H8ATR)  (7) 

Here,  hH2AiR  is  the  molar  flow  rate  of  H2  contained  in  the  ATR 
outlet  stream,  ‘4’.  l(hCOATR)  is  the  molar  flow  rate  of  H2  that  could 
be  produced  from  CO  by  the  WGS  reaction  contained  in  the  ATR 
outlet  stream,  ‘4’.  4(hCH4  ATR)  is  the  molar  flow  rate  of  H2  that  could 
be  produced  from  1  mol  of  CH4  contained  in  stream  ‘4’  (1  mol  of 
CH4  produces  3  mol  of  H2  and  1  mol  of  CO  that  takes  part  in  the 
WGS  reaction  and  converted  to  H2).  The  similar  statement  applies 
for  the  higher  hydrocarbons  e.g.  hC2H6  ATR  and  hCsHg  ATR. 

The  amount  of  H2  and  02  consumed  in  the  anode  stack  is 
calculated  using  the  known  fuel  and  air  utilization  factors, 
respectively  in  an  Aspen  Plus  calculator  block. 

,,  ^H2  consumed  ,  . 

U/  "  * — +u -  (y) 

1  nH2eq 

iiOj, consumed  =  0.5  •  nH2  consumecj  =  0.5  •  Uf  •  flH2  eq  (9) 


Anode  half  -  cell  reaction  :  H2  +  02  ->H20  +  2e 


(10) 


Overall  reaction  :  H2  + 1  /2  02  ->■  H20  (6) 

As  Aspen  Plus  cannot  model  the  electrochemical  reactions 
shown  in  Eqs.  (4)  and  (5),  the  overall  reaction  (Eq.  (6))  is  considered 
instead.  The  stream  ‘C-AIR’,  is  compressed  in  the  compressor  block 


The  calculator  block  calculates  the  equivalent  H2  flow  rate  from 
Eq.  (7).  Once  the  equivalent  H2  flow  rate  is  known  the  consumed  H2 
flow  rate  can  be  calculated  from  Eq.  (8).  The  oxygen  flow  rate  is 
related  to  the  consumed  H2  flow  rate  via  the  overall  electro¬ 
chemical  reaction  (Eq.  (6))  and  can  be  calculated  by  using  Eq.  (9). 
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Aspen  Plus  calculator  block  varies  the  air  flow  rate  to  obtain  the 
desired  oxygen  split  fraction  which  is  also  equivalent  to  the  known 
air  utilization  factor,  Ua.  The  temperature  of  the  depleted  air,  stream 
‘16’  is  raised  by  using  heat  liberated  from  the  electrochemical 
reactions  at  the  heater  module  named  ‘HEATR-4’.  The  stream,  ‘10’, 
still  contains  substantial  amount  of  heat  that  is  calculated  adding 
another  heater  module,  ‘HEATR-3’  by  setting  up  the  temperature 
(450  °C)  of  the  final  exhaust  stream,  ‘11’. 


The  polarization  losses  are  modeled  using  Eqs.  (17)— (26). 
Ohmic  loss  [15] 


a  jpa(AirDm)2 
1  ohm  8  iff 


2.3.  Electrochemical  modeling 

The  cell  voltage  is  calculated  by  applying  the  Nemst  equation 
and  then  subtracting  the  ohmic,  activation,  and  concentration 
losses.  Average  partial  pressures  have  been  used  in  the  Nemst 
equation,  given  in  Eq.  (11). 

l"> 

Ohmic  losses  are  calculated  by  using  equations  proposed  by 
Song  et  al.  [15]  for  tubular  SOFC.  The  activation  losses  are  calculated 
by  assuming  linear  current  potential  equations  for  tubular  cells 
[16].  Semi-empirical  formulas  proposed  by  Achenbach  [17]  are  used 
to  calculate  the  exchange  current  densities  for  activation  voltage 
loss.  Fields  law  of  diffusion  is  used  to  calculate  the  concentration 
voltage  loss  of  anode  and  cathode  for  tubular  cells  considering  both 
ordinary  and  Knudsen  diffusion  phenomena  in  the  form  of  effective 
diffusion  coefficient  and  discussed  briefly  below.  Knudsen  diffusion 
coefficient  for  straight  and  round  pores  [6,16]  is  given  in  Eq.  (12). 


DKi  =  97.0rp . 


(13) 


1  x  1°  ^(l/M.  +  l/M,)’72 

°ik  ~  77:  7772 

''(*»  +*») 

The  effective  diffusion  coefficient  can  be  calculated  from  Eq.  (14). 
Di(eff)  £\DK,i  Dik)  1  ■ 


-A\A  +  2(1  -A-  B)]  (18) 


Vlhm  =J-Pele 


Chm  =  I' Pint' (*■  Dm)'  ' 


Activation  Loss  [16] 


(12) 


Fuller  et  al.  [18]  method  is  considered  to  simulate  the  binary 
diffusion  coefficient 


Ractc  RgT0'Kc\P°2/P°)  eXPL- 


Concentration  polarization  [6,16] 


(19) 

(20) 


(21) 

(22) 


1  -/VVjA 

(  l°  )  j] 

Vs  2  F) 

\Da(eff)  yn2  '  PSOFC  J 

l  +  ^+e.F) 

(  l° 

\Da(e ff) ' Th20  '  PsOFC  J 

The  anode  and  cathode  effective  diffusivities  are  [16]  given  in  Eqs. 
(15)  and  (16). 

Da(eff)  =  Gw)  'DH2(eff)  +  (psofc)  'DH2°(eff)  (15) 


So,  =7 


UK,  02(eff) 


DK,  Oa  (eff)  +  D02-N2  (eff) 

Finally,  the  actual  operating  voltage  of  a  cell, 

V  =  VN-  (rjohm  +  t]3Ct  +  t]con)  (26) 

An  Aspen  Plus  design  specification  block  determines  the  input 
fuel  flow  rate  required  to  achieve  the  desired  DC  power  (Pei,DC  =  VI) 


Table  1 

Composition  of  JP-8  surrogate  fuel  [19], 
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Iso-octane  C8His 

n-Decane  CI0H22 

n-Dodecane  C12H26 

n-Tetradecane  C14H30 

n-Hexadecane  C16H34 

Methylcyclohexane  C7H14 

Cyclooctane  C8H]6 

n-Butylbenzene  Ci0Hi4 

1 ,2,3,4-tetramethylbenzene  C10H14 
m-Xylene  CsHio 

Teralin  C10H12 

1-Methylnaphthalene  CnH10 


formula  Weight  (%)  LHVCkJmor1) 
X66  5064.31 

16.08  6293.01 

22.54  7512.24 

16.87  8731.15 

12.22  9949.55 

3.51  4256.33 

4.54  4912.24 

4.72  5563.34 

4.28  5505.35 

3.95  4330.97 

4.14  5356.48 

3.49  5594.73 


while  calculating  the  voltage  loss  simultaneously.  For  current 
density  sensitivity  analysis  a  calculator  block  determines  the  input 
fuel  flow  rate  based  on  required  current,  I  using  Eq.  (27) 


u  60-/  r-nt 

nH2,eq  =  2P^  (27) 


The  gross  and  net  efficiencies  of  the  system  can  be  calculated  from 
Eqs.  (28)  and  (29),  respectively. 


Pel.AC 

^el’gr0SS  "fuelLHVfuel 

(28) 

Pel, AC  —  Pcomp 

Mnet  ~  tifuel '  LHVfue, 

(29) 

The  composition  of  JP-8  surrogate  fuel  used  in  this  study  is 
shown  in  Table  1. 


Table  2 

Model  input  parameters. 

Geometry  [6,22,23] 

Cell  Length/diameter  (m) 

Anode  thickness,  la  (m) 

Cathode  thickness,  lc  (m) 

Electrolyte  thickness,  le  (m) 
Interconnect  thickness,  ljnt  (m) 
Interconnect  width,  ijnt  (m) 

Material  Properties  [22,24] 

Anode  resistivity,  pA  (n  m) 

Cathode  resistivity,  pc  (£2  m) 
Electrolyte  resistivity,  pe  (O  m) 
Interconnect  resistivity,  pint  (£1  m) 
Ohmic  Loss  [15] 

A/B 

Activation  Loss  [17,25] 
Pre-exponential  factor,  KJKC  (A  nr2) 
Slope,  m 

Activation  energy,  Ea/Ec  (J  mol-1) 
Concentration  Loss  [6,26] 

Electrode  pore  radius, rp  (m) 
Electrode  porosity,  e/tortuosity,  f 


1.5/0.022 

0.0001 

0.0022 


0.000085 

0.009 

2.98  x  10  5exp(  - 1 392/T„p) 
8.114  x  10-5exp(600/Top) 
2.94  x  1 0~5exp( 1 0,350/Top) 
0.025 


2.13  x  10s/1.49  x  1010 
0.25 

100,000/160,000 

5  x  10~7 
0.5/5.9 


it  is  cited  in  numerous  studies  and  the  material  and  geometric 
properties  are  available  to  conduct  the  modeling  study  as  well  as 
performing  validation.  Tables  2  and  3  show  the  model  input 
parameters  and  input  conditions,  respectively.  Table  4  provides 
a  comparison  of  the  present  study  with  operation  and  previous 
simulation  data.  The  present  simulation  results  matched  very 
closely  with  the  simulation  results  published  by  Doherty  et  al.  [6] 
but  differs  from  the  experimental  plant  data  due  to  the  simplified 
assumptions  in  this  study.  Nevertheless,  the  error  associated  with 
the  current  density  (j)  and  gross  AC  efficiencies  are  1.14%  and  3.34%, 
respectively,  compared  to  the  operation  conditions. 


3.  Model  validation 


4.  Model  results 


The  ATR  module  was  individually  validated  against  experi¬ 
mental  data  published  by  Liu  et  al.  [20]  of  a  lcW  scale  catalytic  ATR 
using  hexadecane  and  dodecane  as  a  diesel  surrogate  fuel.  The 
simulation  result  matched  more  closely  with  the  experimental 
results  at  higher  O2/C  ratios  shown  in  Fig.  2.  Less  fuel  reforming  is 
likely  to  occur  at  lower  O2/C  ratio,  resulting  in  insufficient  heat  to 
promote  the  steam  reforming  reactions,  which  may  be  the  reason 
for  this  difference  at  lower  O2/C  ratios. 

To  validate  the  tubular  SOFC  module  in  this  study  a  Siemens 
Power  Generation  Inc  (SPGI)  100  kW  natural  gas  fueled  cogenera¬ 
tion  system  was  selected.  The  100  kW  SPGI  system  is  the  world’s 
longest  operational  SOFC  generator  with  more  than  36,900  h 
(updated  on  June  2007)  of  operation  [21],  This  model  is  selected  as 


Fig.  2.  Comparison  of  experimental  (symbols)  reformate  concentrations  [20]  with 
Aspen  Plus  simulation  (lines)  for  dodecane  at  H20/C  =  2.0  with  varying  02/C. 


4.1.  Setting  the  ATR  inlet  air  temperature 

The  ATR  air  inlet  temperature  was  varied  from  350  °C  to  550  °C 
over  a  H2O/C  ratio  range  of  0.1— 1.0,  to  determine  the  optimal  inlet 
air  temperature.  These  simulations  were  conducted  at  three 
reformer  operating  temperatures  e.g.  700  °C,  750  °C  and  800  °C. 
Figs.  3  and  4  show  the  effect  of  ATR  air  inlet  temperature  at 
a  representative  ATR  operating  temperature  of  750  °C.  With 
increasing  reformer  air  inlet  temperature  the  O2/C  ratio  decreased 


Table  3 

Model  input  conditions  for  tubular  SOFC  based  system  validation  [5,6,27], 


Fuel  composition,  (X  vol.) 


Operating  pressure  (PsoFc)/Ejector 
pressure  ratio 

Cell  operating  temperature,  Top  (°C) 

Inlet  fuel  temperature,  (°C) 

Inlet  air  temperature,  (°C) 

H2O/C 

Air  utilization  factor,  Ua 
Fuel  utilization  factor.  Up 
Ejector  fresh  fuel  pressure  ratio 
Temperature  difference  between  the  outlet 
of  cold  and  hot  stream  of  HEATX,  (°C) 
Power  Output  DC,  Pei,rx:  (kW) 

DC  to  AC  inverter  efficiency,  (%) 

Pressure  drops  inside  the  SOFC,  (bar) 

SOFC  thermal  loss,  (%) 


CH4  81.3,  C2H6  2.9,  C3H8  0.4, 
C,H10  0.2,  N2  14.3,  C02  0.9 
109,431/3 

96.0768  (1152  cells) 

910 

200 

630 

1.8 

0.19 

0.85 

3.0 

10 

120 

92 
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Table  4 

Model  results  compared  to  literature. 


Voltage,  (V) 

Current  density,  (mA  cnr2) 
Pre-reforming  temperature,  (°C) 
Pre-reformer  CH4  conversion,  (%) 
Cathode  inlet  temperature,  (°C) 
Combustion  products  temperature,  (°C) 
Stack  exhaust  temperature,  (°C) 

Anode  inlet  gas  composition,  (mole  %) 

Anode  outlet  gas  composition,  (mole  %) 

Stack  outlet  gas  composition,  (mole  %) 
Gross  AC  efficiency  (LHV),  (%) 

Net  AC  efficiency  (LHV),  (%) 
nr  =  not  reported. 


Operation  data  [5,23,27,28] 


550 


H2  14.0,  CO  5.0,  H20  48.0, 

C02  28.0,  N2  5.0 

N2  77.0,  02  16.0,  H20  5.0,  C02  2.0 
50 


Simulation  data  [6] 


182.86 

535.1 

25 

823.7 
1012.3 

833.7 

H2  26.9,  CO  5.6,  CH4  10.4, 

H20  27.8,  C02  23.1, N2  6.2 
H2  11.6,  CO  7.4,  H20  50.9, 

C02  24.9, N2  5.1 

H20  4.5,  C02  2.3, 02  15.9,  N2  77.3 

51.28 

49.15 


1016 


H2  26.9,  CO  5.6,  CH4  10.4,  H20  27.7, 

C02  23.1, N2  6.2 

H2  11.6,  CO  7.4,  H20  50.9,  C02  24.9, N2  5.1 


H20  4.7,  C02  2.4, 02  15.6,  N2  77.2 

51.67 

48.42 


in  the  H20/C  ratio  range  of  0.1-1.0.This  was  a  result  of  the  ATR  air 
conveying  more  thermal  energy  into  the  reformer  subsequently 
requiring  less  oxygen  for  adiabatic  operation,  and  thus  lowering 
fuel  dilution  by  N2.  At  higher  ATR  air  inlet  temperature  less  oxygen 
is  required  by  the  reformer,  reducing  direct  oxidation  of  H2  and  CO 
(results  not  reported).  Based  upon  these  observations  the  reformer 
inlet  air  temperature  of  550  °C  is  chosen  to  reduce  the  number  of 
simulations  associated  with  the  system. 


4.2.  Effect  of  H2O/C  on  ATR  product  gas  composition 

Reforming  higher  hydrocarbons  requires  preventing  of  coke 
formation  to  avoid  the  blockage  of  catalyst  pores  and  reaction 
active  sites  that  may  eventually  make  the  reformer  inoperative.  The 
anode  recycle  also  contains  C02,  CO  and  H2,  as  modeled  by  Eq. 
(30)— (33). 

C  +  C02<-  2CO,  AH29815  =  +172.422  kj/mol, 

AG  <  0  at  T>700  °C  t30) 

C  +  H20~C0  +  H2,  AH29815  w*131.285kJ/mol, 

AG  <  0  at  T>674  °C 

CO  +  3H2 «CH4  +  H20,  A H298  15  =  — 205.8kJ/mol, 

AG  <  0  at  T>620  °C  [  j 


C0  +  H20«C02  +H2i  A H29815  =  -41.138kJ/mol, 

AG  <  0  atT>816.5  °C  L>J’ 

Fig.  5  shows  the  effect  of  H20/C  ratio  on  the  reformer  outlet  gas 
composition  at  a  representative  reformer  temperature  of  700  °C. 
The  02/C  ratio  quickly  reached  a  maximum  value  at  H20/C  ratio  of 
0.4,  followed  by  a  steady  decline  with  higher  H20/C  ratios.  At  low 
H20/C  ratios,  the  reformer  required  more  oxygen  to  provide  the 
energy  necessary  for  adiabatic  operation.  However,  as  the  H20/C 
ratio  increased,  anode  recycle  also  increased,  reducing  the  reformer 
oxygen  requirement. 

Increased  H20/C  ratio  requires  higher  anode  recycle  fraction 
bringing  more  steam  and  C02  into  the  reformer  taking  part  in 
the  reverse  Boudouard  reaction  (Eq.  (30))  and  carbon  steam 
reforming  reaction  (Eq.  (31))  which  are  thermodynamically 
favorable  at  this  temperature  range.  As  a  result,  solid  carbon  (in 
the  form  of  graphite)  formation,  C(gr)  decreases  rapidly  with 
H20/C  ratio  shown  in  Fig.  5.  The  reverse  methanation  reaction 
(Eq.  (32))  also  becomes  thermodynamically  favorable  at 
temperatures  greater  than  620  °C  and  helped  suppress  methane 
formation  with  H20/C  ratio.  CO/C  ratio  reached  a  peak  value  of 
1.046  with  increasing  H20/C  due  to  the  favorability  of  reverse 
Boudouard  reaction,  carbon  steam  reforming  reaction,  and 
reverse  methanation  reactions  in  the  coke  formation  region. 
Solid  carbon  may  also  take  part  in  the  direct  oxidation  reaction 
in  the  presence  of  oxygen  and  thus  can  produce  more  CO.  The 
H20/C  ratio  at  which  solid  carbon  formation  becomes 


T.  Tanim  et  al.  /  Journal  of  Power  Sources  221  (2013)  387—396 


393 


Fig.  5.  Mole  flow/C  at  the  ATR  outlet  stream  for  tubular  SOFC  system  at  Tre f  =  700  °G 


thermodynamically  unfavorable  is  the  minimum  H2O/C  ratio 
above  which  the  system  should  operate. 

H2  production  increases  with  H2O/C  ratio  due  to  the  occurance 
of  more  reforming  reactions  with  the  JP-8  fuel  and  WGS  reaction 
(Eq.  (33)).  But  the  rate  of  H2  production  was  not  strongly  dependent 
upon  the  WGS  reaction  due  to  its  weak  thermodynamic  favorability 
at  this  elevated  temperature.  This  resulted  in  marginal  CO  decrease 
in  the  coke-free  region.  Some  H2  may  have  also  been  oxidized  due 
to  the  presence  of  oxygen  in  the  ATR,  producing  steam. 

4.3.  Effect  of  ATR  temperature 

ATR  performance  was  quantified  at  four  operating  tempera¬ 
tures:  700  °C,  750  °C,  800  °C,  and  850  °C.  Solid  carbon  formation 
reduces  substantially  with  increasing  ATR  operating  temperature 
due  to  the  high  thermal  decomposition  of  the  JP-8  fuel  and  the 
high  probability  of  the  reverse  Boudouard  and  carbon  steam 
reforming  reactions  shown  in  Fig.  6.  Favorability  of  reverse 
methanation  reaction  in  the  presence  of  higher  H2O/C  ratio  also 
prevents  direct  thermal  cracking  of  methane  thus  reducing  coke 
deposition.  The  minimum  H2O/C  ratio  to  avoid  any  coke  formation 
calculated  are  0.4,  0.25,  0.15  and  0.12  for  ATR  operating  temper¬ 
atures  of  700  °C,  750  °C,  800  °C,  and  850  °C,  respectively. 
Production  of  H2  decreases  with  temperature  while  production  of 
CO  increases  shown  in  Figs.  7  and  8.  At  higher  temperatures,  the 
WGS  reaction  weakens  and  becomes  thermodynamically  unfa¬ 
vorable  at  816  °C.  Also,  the  anode  recycle  contains  substantial 
amount  of  CO2  which  promotes  reverse  WGS  reaction  at  high 
temperatures  to  produce  additional  CO.  Additionally,  direct 
oxidation  of  H2  to  steam  at  high  temperatures  might  have  been 
another  cause  for  lower  H2  production. 


h2o/c 


Fig.  7.  H2/C  ratio  at  different  ATR  operating  temperatures. 


h2o/c 

Fig.  8.  CO/C  at  different  ATR  operating  temperatures. 


4.4.  Performance  of  the  SOFC  system 

The  performance  of  the  modeled  tubular  SOFC  based  system 
maintained  at  a  constant  temperature  of  910  °C  was  investigated  at 
Ua  —  0.167  and  Uf=  0.85  in  the  coke-free  region.  Results  presented 
in  Fig.  9  shows  that  cell  voltage  decreased  with  H2O/C  ratio, 
possibly  due  to  the  decrease  of  H2  and  CO  concentrations  at  the 
anode  outlet  stream.  Higher  H20/C  ratios  required  a  higher  anode 
recycle  split  fraction  (volume  basis)  to  meet  the  steam  and  energy 


394 


T.  Tanim  et  at  /  Journal  of  Power  Sources  221  (2013)  387-396 


h2o/c 

Fig.  10.  Cell  voltages  at  different  ATR  operating  temperatures. 


requirements  of  the  reformer.  Additional  steam  from  the  anode 
recycle  led  to  higher  H2  production  through  greater  shifting,  but  the 
overall  concentrations  of  H2  and  CO  gradually  decreased  due  to  the 
high  volumetric  flow  rate  of  recycle  stream,  causing  a  decay  in 
Nernst  cell  voltage.  Current  density  increased  proportionately  with 
lower  voltage  and  with  H20/C  ratio,  while  the  ohmic,  activation  and 
concentration  losses  increased  gradually  as  well.  The  cathode  air 
flow  rate  increased  to  provide  more  oxygen  in  order  to  support  the 
electrochemical  reactions  at  the  SOFC  stack.  The  net  AC  efficiency 
decreased  due  to  the  higher  parasitic  power  required  by  the 
cathode  air  compressor  to  provide  this  additional  air  requirement. 

The  H2  and  CO  concentrations  at  the  ATR  outlet  stream/SOFC 
inlet  stream  were  lower  at  higher  ATR  temperatures  as  discussed  in 
the  previous  sections,  thus  reducing  cell  voltage  as  shown  in  Fig.  10. 
Increasing  ATR  temperature  raised  the  average  temperature  of  the 
SOFC  streams  thus  lowering  Gibb’s  energy  and  the  resulting  Nernst 
cell  voltage.  Nevertheless,  the  system  showed  the  highest  cell 
voltage  of  671.7  mV  at  700  °C  ATR  operating  temperature  and  H20/ 
C  =  0.4  and  the  lowest  cell  voltage  of  567.1  mV  at  850  °C  ATR 
operating  temperature  and  H20/C  =  1.0.  As  shown  in  Fig.  11,  net  AC 
efficiency  decreased  with  ATR  operation  temperature  due  to 
parasitic  loss  associated  with  higher  ATR  and  cathode  inlet  air  flow 
rates  and  higher  fuel  flow  rate  (not  presented)  required  to  achieve 
the  specified  power  output. 

Fig.  12  presents  the  total  efficiency  (summation  of  net  AC  effi¬ 
ciency  and  thermal  efficiency)  of  the  system  at  different  ATR 


Fig.  12.  Total  efficiency  at  different  ATR  operating  temperatures. 


Fig.  13.  Effect  of  Uf  on  voltage,  efficiency  and  current  density  at  Tref  =  700  °C  and  H20/ 


HJO/C 

Fig.  11.  Variation  of  net  AC  efficiency  at  different  ATR  operating  temperatures. 


Fig.  14.  Effect  Uf  on  post  combustion  temperature  and  JP-8  flow  rate  at  Tref  =  700  °C 
and  H20/C  =  0.5. 


operating  temperatures.  The  thermal  efficiency  has  been  calculated 
keeping  the  temperature  of  the  exhaust  (stream-11 )  at  450  °C.  It  is 
assumed  that  the  remaining  exhaust  heat  can  be  used  to  vaporize 
the  fresh  fuel.  Total  efficiency  was  higher  at  lower  H20/C  ratio  due 
to  higher  thermal  efficiency  attributed  by  availability  of  more 
unburned  H2  and  CO  at  the  anode  exhaust  stream. 

4.5.  Sensitivity  analysis  of  the  5  kWe  system 

Sensitivity  analyses  were  performed  for  the  5  kWe  system  at  the 
highest  and  lowest  ATR  temperatures  to  show  the  influence  of  fuel 
utilization  factor,  Uf  and  the  current  density,  j  in  the  coke-free 
regions.  700  °C  and  850  °C  ATR  operating  temperatures  were 
selected  for  analysis  as  the  system  demonstrated  highest  and 
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Fig.  15.  Effect  of  current  density  on  power  density,  efficiency  and  voltage  at 
Tref  =  700  "C  and  H20/C  =  0.5. 


lowest  net  AC  efficiency,  respectively  at  these  temperatures.  H2O/C 
ratios  of  0.5  and  0.3  were  maintained  corresponding  to  the  ATR 
operation  temperatures  of  700  °C  and  850  °C,  respectively.  These 
H2O/C  ratios  are  slightly  greater  than  the  minimum  H2O/C  ratio 
required  to  prevent  coke  formation  in  the  maximum  efficiency 
zone.  The  air  utilization  factor,  Ua  was  kept  constant  at  0.167  by 
varying  the  air  flow  rate. 


Fig.  17.  Effect  of  current  density  on  stack  AC  power,  air  flow  rate  and  fuel  flow  rate  of 
the  system  at  Tref  =  700  °C  and  H20/C  =  0.5. 


45.1.  Fuel  utilization  factor,  Uf 

Data  presented  in  Fig.  13  indicates  that  cell  voltage  decreased 
slightly  due  to  reduced  H2  content  at  the  anode  outlet  stream 
and  more  anode  polarization  loss.  Current  density  increased 
with  Uf  as  more  H2  was  consumed  (1  =  2F  x  nH2, consumed)  at  the 
cells  with  higher  Uf.  The  ohmic,  activation  and  concentration 
losses  are  proportional  to  the  current  density  causing  a  higher 
voltage  loss  with  higher  Uf.  At  low  utilization  factors  more 
unburned  H2  and  CO  go  into  the  post  combustor  and  release 
more  chemical  energy  in  the  form  of  heat,  increasing  the 
thermal  efficiency  of  the  system.  This  causes  a  large  temperature 
increase  at  the  post  combustor  outlet,  shown  in  Fig.  14.  The  high 
post  combustion  temperature  increased  the  cathode  inlet  air 
temperature,  which  decreased  the  available  Gibb's  free  energy 
and  thus  reduced  Nernst  voltage.  The  net  AC  efficiency  increased 
as  it  is  proportional  to  both  Uf  and  cell  voltage.  Increasing  the 
fuel  utilization  factor  implies  more  H2  would  be  consumed  by 
the  SOFC  stack  and  therefore  converting  more  chemical  energy 
into  electrical  energy,  requiring  less  fresh  JP-8  fuel,  as  shown  in 
Fig.  14. 

4.5.2.  Current  density,  j 

Fig.  15  demonstrates  the  effect  of  current  density  on  cell 
voltage,  power  density,  gross  and  net  AC  efficiencies.  It  is  observed 
that  both  voltage  and  efficiency  decreased  with  current  density 
while  power  density  increased.  As  shown  in  Fig.  16,  cell  voltage 
decreased  significantly  with  current  density  because  of  higher 
voltage  losses  that  are  directly  proportional  to  j.  The  concentration 
loss  increased  sharply  at  higher  current  densities  due  to  the 
inability  of  gaseous  species  to  diffuse  rapidly  through  the  porous 
electrodes  at  these  high  current  densities.  Ohmic  loss  was  domi¬ 
nant  in  the  tubular  SOFC  based  system  due  to  the  long  current 
path  of  electron  travel  around  the  tubular  cells.  The  systems 
showed  consistent  results  at  other  ATR  operating  temperatures  as 

The  net  AC  efficiency  of  the  system  decreased  due  to  the 
parasitic  loss  associated  with  the  compression  of  cathode  and 
oxidant  air  streams  whose  flow  increased  with  current  density 
shown  in  Fig.  17.  In  this  analysis,  fuel  flow  rate  increased  at  higher 
current  density,  which  also  lowered  the  system  efficiency.  Power 
density  increased  to  a  peak  value  of  137.8  mW  cm-2  and 
134.7  mW  cm“2  at  a  current  density  of  300  mA  cm’2  for  ATR 
operating  temperatures  of  700  °C  and  850  °C,  respectively.  In 
terms  of  efficiency,  the  data  indicates  that  it  was  beneficial  to 
operate  the  cells  at  lower  current  densities  to  achieve  higher 
efficiency,  thus  reducing  operating  cost  and  lowering  power 
production.  So  to  obtain  the  desired  power,  an  actual  system  will 
require  more  cells  making  the  system  relatively  larger  and  more 
capital  cost  intensive.  On  the  other  hand,  operating  an  actual 
system  at  higher  current  densities  would  require  less  number  of 
cells,  making  the  system  less  expensive  but  efficiency  will  be 
lower  causing  higher  operating  cost.  A  trade-off  between  oper¬ 
ating  cost  and  capital  cost  based  on  voltage,  efficiency  and  power 
needs  to  be  conducted  to  set  the  system  at  an  optimum  operating 
point.  Table  5  shows  recommended  operating  conditions  based  on 
this  study. 


Recommended  operating  current  density,  j  for  the  systems. 

rref(°C)  H20/C  j(mAcm-2)  Cell  voltage  (V)  tjei_net(%)  Pac(kW)  Power  density  (mW  cm-2) 

700  05  175-225  0.68-0.6  40.6-36.2  4.85-5.58  109.8-126.3 

850  0.3  0.67-0.6  36.6-32.6  4.78-5.48  108.1-124.1 
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5.  Conclusion 

A  desulfurized  JP-8  fuel  based  5  kWe  SOFC  power  unit  has  been 
modeled  in  this  study  making  reference  to  the  performance  of  tubular 
cells  developed  by  SPGI.  Results  indicated  that  characterization  of  the 
ATR  is  a  critical  step  in  system  evaluation,  as  SOFC  stack  performance 
will  be  largely  dependent  upon  its  operation.  Further,  H2O/C  ratio  is 
another  key  system  operating  parameter,  because  it  governs  coke 
formation  and  H2  production.  This  study  found  a  minimum  H2O/C 
ratio  of  0.4,  0.25,  0.15,  and  0.12  needs  to  be  maintained  for  coke-free 
reforming  at  ATR  operating  temperatures  of 700  °C,  750  °C,  800  °Cand 
850  °C,  respectively.  But  increasing  H20/C  ratio  causes  fuel  dilution 
which  reduces  the  cell  voltage  and  net  system  efficiency.  Superior 
system  performance  has  been  observed  at  lower  ATR  operating 
temperature  and  lower  H2O/C  ratio  in  the  coke-free  region.  The  system 
demonstrated  a  maximum  efficiency  of  39.5%  at  700 0  C  ATR  operating 
temperature  and  H20/C  =  0.4  and  minimum  of  35.9%  at  850  °C  ATR 
operating  temperature  and  H20/C  =  0.12  in  the  lower  H20/C  zone.  In 
the  higher  H2O/C  ratio  zone  the  maximum  net  AC  efficiency  of  the 
system  has  always  been  observed  at  H2O/C  =  1.0  which  was  39.5%  and 
32.6%  at  700 0  C  and  850 0  C  ATR  temperatures,  respectively.  The  system 
demonstrated  a  maximum  total  efficiency  of  72.4%  at  700  °C  ATR 
operating  temperature  with  H2O/C  ratio  of  0.4  and  minimum  of  65.6% 
at  850  °C  ATR  operating  temperature  with  H20/C  ratio  of  1.0. 

Fuel  utilization  is  another  key  parameter  identified  through  this 
study.  The  parametric  analyses  showed  the  cell  voltage  marginally 
decreases  with  17/ but  an  increase  in  net  AC  efficiency  was  observed. 
Total  efficiency  decreases  with  Uf  primarily  due  to  reduction  of 
thermal  efficiency.  Greater  power  density  can  be  achieved  with  higher 
current  density  operation  but  at  the  expense  of  cell  voltage  and  net  AC 
efficiency.  The  system  showed  a  maximum  efficiency  of  40.6—36.2% 
and  36.6-32.6%  at  the  current  density  range  of  175-225  mA  cm2 
and  ATR  operating  temperatures  of  700  °C  and  850  °C,  respectively, 
corresponding  to  H2O/C  ratio  of  0.5  and  03,  respectively. 
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